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Abstract. 

Using the explicit numerical solution of the axially-symmetric Gross-Pitaevskii 
equation we study the dynamics of interaction among vortex solitons in a rotating 
matter-wave bright soliton train in a radially trapped and axially free Bose-Einstein 
condensate to understand certain features of the experiment by Strecker et al. [2002 
Nature 417 150]. In a soliton train, solitons of opposite phase (phase 5 = tt) repel 
and stay apart without changing shape; solitons with 6 = attract, interact and 
coalesce, but eventually come out; solitons with a general 6 usually repel but interact 
inelastically by exchanging matter. We study and suggest future experiments with 
vortex solitons. 
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1. Introduction 

Solitons are solutions of wave equation where localization is obtained due to a nonlinear 
interaction. Solitons have been noted in optics jTj, high-energy physics and water 
waves |2j, and more recently in Bose- Einstein condensates (BEC) [SlIllE]- The bright 
solitons of BEC represent local maxima jHElini, whereas dark solitons represent local 
minima [SlIZ]- In addition to the observation of an isolated bright soliton in an expulsive 
potential [Sj , a number of bright solitons constituting a soliton train was also observed 
in an experiment by Strecker et al. [4 , where they suddenly turned a repulsive BEC 
of ^Li atoms attractive. Consequently, the BEC collapsed, exploded and generated a 
soliton train. Similar collapse and explosion in ^^Rb BEC have been studied before from 
a different point of view 0. These experiments were performed by manipulating the 
background magnetic field near a Feshbach resonance It was found |4j that solitons 
in such a train usually stay apart. Also, often a soliton was found to be missing from 
a train |3]. Hasegawa considered the generation of a train of optical solitons by an 
induced modulational instability. There have been theoretical attempts ^11V12"'1H"'14 IITK] 
to simulate essentials of these experiments jHE] on bright solitons in BEC. Al Khawaja 
et al. jTT] used a variational approach to describe the experiment by Strecker et al. |lj 
whereas Salasnich et al. [T2j and Leung et al. [TH] used an effective one-dimensional 
model for the same purpose. In a more recent work Salasnich et al. [14J considered a 
three-dimensional mean-field model to study the production of and interaction between 
bright solitons to account for different aspects of the experiment by Strecker et al. |3]. 
Carr et al. [15^ used an approximate approach to describe the experiment of Khaykovich 
et al. Elyutin et al. and Shchesnovich et al. studied solitons in one dimension. 

We use the explicit numerical solution of the axially- symmetric mean-field Gross- 
Pitaevskii (CP) equation JTj to study the dynamics of bright solitons in a soliton 
train It is found that in an axially-symmetric configuration with no axial trap, the 
GP equation permits solution in the form of a radially confined soliton train. We also 
give an explanation of missing solitons in the experiment of Strecker et al. 

The present approach is also extended to the study of vortex solitons with an 
angular momentum of h per atom in the axial direction. Vortex solitons are rotating 
solitons of an attractive condensate. Due to experimental observation JH] of a vortex 
state in a rotating BEC, the experimental generation of vortex solitons seems possible. 
Attractive BEC's may not form vortices in a thermodynamically stable state. However, 
such vortices may be created via a Feshbach resonance jH]- Due to the conservation 
of angular momentum, a vortex soliton train could be generated by suddenly changing 
the inter-atomic interaction in an axially-symmetric rotating vortex condensate from 
repulsive to attractive near a Feshbach resonance in the same fashion as in the 
experiment by Strecker et al. jl] for a non-rotating BEC. Alternatively, a single vortex 
soliton could be prepared and studied in the laboratory by forming a vortex in a 
small repulsive condensate and then making the interaction attractive via a Feshbach 
resonance and subsequently reducing the axial trap slowly. Already, there have been 
theoretical considerations for these vortex states JHH^Hj- It would be worthwhile to 
explore these possibilities experimentally. 

In particular we study in some detail the interaction between two bright vortex 
solitons for different relative phases between them. Usually, the interaction is found to 
be repulsive and inelastic with exchange of particles. The interaction turns attractive 
for small values of relative phase. In one dimension the interaction is usually elastic 
without exchange of particles |2j. 

In section 2 we present the mean-field model of solitons that we use in the numerical 
analysis. In section 3 we present the results of our numerical study. Finally, in section 
4 we present the conclusions. 
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2. Nonlinear Mean-free Model for Solitons 



In the following one- dimensional nonlinear free Schrodinger equation in dimensionless 
units 

|vl>(y,t)|2 vl/(^,t) = (2.1) 



d 



dt dy'^ J 

solitons are bound states due to the attractive nonlinear interaction with wave function 



at time t and position '^{y,t) = y/2|n[exp(—zfit)sech(?/A/|f^|), with f2 the energy 
Equation (1) can sustain any number of such solitons at different positions. Although, 
the one-soliton solution of (1) may remain stationary at a fixed position, the many 
solitons of a multiple-soliton solution of this equation generally move around because of 
the interaction among them 

The time-dependent Bose- Einstein condensate wave function \E'(r;r) at position r 
and time r is described by the following mean-field nonlinear GP equation J7] 
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where m is the mass and N the number of atoms in the condensate, g = Anh'^a/m the 
strength of inter-atomic interaction, with a the atomic scattering length. For an axial 
trap V^(r) = |ma;^(r^ -|- A^z^) where a; is the angular frequency in the radial direction 
r and Xu that in the axial direction z, with A the aspect ratio. The normalization 
condition is J (ir|\E'(r; r)p = 1. 

In a quantized vortex state with each atom having angular momentum Lh 
along the z axis, \E'(r, r) = ilj{r, z,t) exp{iL6) where 6 is the azimuthal angle. Now 
transforming to dimensionless variables x = y/2r/l, y = y/2z/l, t = tu, I = ^/fij {moo), 
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(2.3) 



where non-linearity n = Na/l. For solitonic states n is negative. In terms of the 
one-dimensional probability P(?/, t) defined by 



P{y,t) = 2n dx\^{x,y,t)\yx, 



(2.4) 

the normalization of the wave function is given by dyP{y, t) = 1 

We solve the GP equation ()2.3|1 numerically using a variation of the split-step time- 
iteration method using the Crank-Nicholson discretization scheme described recently 
j23j . Typical space and time steps for discretization are 0.1 and 0.001. The variation 
of the standard approach is required for A = 0. The time iteration is started with the 
known harmonic oscillator solution for a small A = Aq ~ 0.05 and nonlinearity n = 0: 
ip{x,y) = [A/{22^+37r3(L!)2}]i/4 ^i+L^-ix^+Xy^)/4 ^^^Yi energy (1 + L + A/2) pOI. The 
desired value of nonlinearity n is then slowly switched on and A slowly switched off 
from A = Ao to zero in the course of time iteration. The solution then corresponds to 
the trapped BEG for A = 0. Then, without changing any parameter, this solution is 
iterated several thousand times so that a solution for A = is obtained independent of 
the initial input A = Aq. 
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3. Numerical Results 

A classic soliton in three-dimensional BEC can be realized for attractive nonlinear 
potential {n < 0) by setting A = [U] in ()2.H|1 . For a fixed L, the BEC is then 
governed by the single parameter n. The absence of a trap in the axial y direction 
will allow a free movement of the solitons in this direction. Consequently, the study of 
soliton interaction will be trap independent. A localized BEC so created should be the 
three-dimensional analogue of the one-dimensional soliton. 

However, for calculational or experimental convenience, in the recent studies some 
weak potential was applied in the axial y direction. In the classic experiment of Strecker 
et al. |3] an optical trap was maintained in the y direction. In the recent theoretical 
investigation by Al Khawaja et al. jTIj on this experiment a weak harmonic trap was 
applied in this direction, whereas Salasnich et al. ^2] employed infinite walls in an 
effective one-dimensional model. In the experiment by Khaykovich et al. [5 and in the 
related theoretical study by Carr et al. [12] an expulsive potential (A^ < 0) was applied 
in y direction. We note that for A^ < 0, only meta-stable and no stable soliton of ()2.3j) 
is possible as the potentials in this equation including the attractive nonlinear term do 
not provide confinement jTHj. For the same reason no bright soliton can be generated 
for n > in ()2.3|) (repulsive condensate). 

Although, under the conditions n < and A = the potentials of ()2.3|) lead 
to confinement, a soliton-type BEC state can be generated only for n greater than a 
critical value (ncr): 'T-cr < n <{]. For n < uqy, the system becomes too attractive and 
collapses and no stable soliton could be generated. The collapse was first confirmed in 
the pioneering experiment by Gerton et al. |24j for a trapped BEC of ^Li. The actual 
value of ncr is a function of the trap parameter A. For the spherically symmetric case 
A = 1, and ncr = —0.575 f^EO]. The value of ncr should slightly change for A = 0. 

We solve the CP equation for solitons with A = and L = and 1 [20] • We 
find numerically that the critical n for collapse of a single soliton is ncr = —0.67 for 
L = and ncr = —2.10 for L = 1. The result for L = is in good agreement with (a) 
the numerical result —0.676 obtained by Gammal et al. using the Crank-Nicholson 
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Figure 2. One-dimensional probability P{y,t) vs. y and t for a train of (a) two 
solitons each with n = —0.2, S — tt, L — 0, and (b) four i = 1 vortex solitons each 
with n = —0.4 and a phase difference S — tt between neighboring solitons. 

method, (b) the numerical result —0.676 obtained by Perez-Garcia et al. ^ using the 
steepest-descent method to minimize the mean-field Hamiltonian [these authors quote 
Q = —SnNa/l = 17 instead of n = Na/l] and (c) the result —2/3 obtained by Salasnich 
et al. [T^ using an approximate analytic one- dimensional model. By solving the three- 
dimensional GP equation Salasnich et al. ^2] found their result to be very accurate. 
However, there is some discrepancy between these results for L = and the value 
'^cr = —0.6268 ± 0.0035 obtained by Carr and Castin fTS] using the imaginary time 
relaxation method. Further independent studies are necessary to resolve the discrepancy. 
For a; = 27r X 800 Hz and final scattering length — 3ao as in the experiment of Strecker et 
al. [T, uqy = —0.67 corresponds to about 6000 ^Li atoms. One can have proportionately 
about three times more atoms in the L = 1 state. 

A L = soliton with n = —0.2 is illustrated in figure 1 (a) where we plot the three- 
dimensional wave function \ip{x,y)/x\ vs. x and y. Next we consider L = 1. In this 
case we calculated the soliton for n = —1 and plot the three-dimensional wave function 
\(p{x,y)/x\ in figure 1 (b). Because of the radial trap the soliton remains confined in 
the radial direction x, although free to move in the axial y direction. In either case the 
single soliton remains stable for more than 400 000 time iterations. 

After having demonstrated the formation of a single soliton we next consider the 
dynamics of two solitons in a soliton train. Two solitons are prepared at positions yi and 
y2 and then superposed with a phase difference 6. Specifically, we consider the following 
superposition of two normalized solitonic waves at y = ±yQ with phase difference o 
between them 

^ix,y) = \(p{x,y + yo)\ + e'^\<f{x,y -yo)\, (3.1) 

with 2yo the initial separation between the solitons in the axial direction. To c onserve 
the total number of atoms we normalize the superposed wave function 1)3.11) of the 
two solitons to 2 as it contains twice as many particles as in a single soliton. In 
this fashion one can also construct the superposition of several solitons with a specific 
phase difference between them. The time evolution of t he s oliton train so formed upon 
superposition is found using the iterative solution of ()2.3|) . We present results for a 
soliton train with two vortex solitons of same nonlinearity with a phase difference 6 of 
7r,37r/4, 7r/2, 7r/4, 7r/8 and between them. Consideration of two equal solitons does 
not lead to a specialization and we shall see that for a general 6 two equal initial solitons 
generally lead to two unequal solitons due to exchange of atoms. 

In one dimension the solitons attract for 6 = and repel for 6 = tt [22lE]- In 
the present three-dimensional case for a general 6 a more complicated motion emerges. 
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In the repulsive 6 = tx case, two solitons stay away from each other. In the attractive 
5 = case, they come close, coalesce, interact and come out. In a one-dimensional 
model of two trapped (A 7^ 0) three-dimensional solitons, it has been shown that these 
solitons repel for (5 = tt [THIT^ . In the following we study the rich dynamics of soliton 
interaction in a bright soliton train in three dimensions. The solitons can easily be set 
into motion by applying a perturbation or a constant (gravitational) force in y direction. 
However, for calculational convenience we chose not to do that and thus we studied the 
interesting relative motion between solitons suppressing their center-of-mass motion. In 
all the cases studied the initial velocity of the solitons is zero. The numerical simulation 
was performed on a lattice 7 > x > and 250 > y > —250. 

We studied the dynamics of interaction of two solitons in view of the experiment 
by Strecker et al. |4j. For simplicity first we prepared two solitons of equal mass 
corresponding to the same nonlinearity n = —0.2 centered at points y = ±15 with 
angular momentum L = and introduced them as the input to the GP equation with 
the initial phase difference S = n. The time evolution of the train of two such solito ns is 
exhibited in figure 2 (a) where we plot the one-dimensional probability P{y, t) of ()2.4j] 
vs. y and t. Because of mutual repulsion for S = tt, the two solitons stay apart and 
move away from each other. In an interval of time 200, the two solitons moved from 
positions |/ = ±15 to ±25, respectively. 

We studied soliton trains with both L = and L = 1 solitons and they lead 
to physically similar results. In the rest of the study we consider only L = 1 vortex 
solitons as they were never studied before. Also, an effective one-dimensional model 
usually employed in other studies does not seem to be applicable to the study of a 
vortex soliton. The rotational degree of freedom responsible for the generation of vortex 
does not exist in a one-dimensional model. We consider four L = 1 vortex solitons each 
with n = —0.4 at positions ±45 and ±15 with phase difference 6 = n between two 
neighboring and hence, repelling, solitons. Obviously, one can accommodate as many 
repelling solitons as one prefers in a train. Because of the repulsion, the solitons stay 
apart and move forward without interaction and maintaining shape. Such a soliton train 
is illustrated in figure 2 (b). In an interval of time 200, the four solitons have moved 
from positions y = ±15 and ±45 to ±17 and ±52, respectively. 

Now we consider the interaction between two vortex solitons in some detail. For 
that we reduced the phase 6 between two L = 1 vortex solitons ai y = ±15 each 
with n = —0.4 from vr to 0. In figures 3 (a), (b), (c), (d), (e) and (f) we plot one- 
dimensional probability P{y, t) of ()2.4|) for two interacting vortex solitons vs. y and t 
for S = IT, 37r/4, 7r/2, 7r/4, 7r/8 and 0, respectively. As S is reduced from vr to 7i/2 the 
two solitons continue to repel and stay apart without much change in shape. However, 
while 6 is reduced from vr the repulsive interaction between the solitons gradually become 
"inelastic" and there is exchange of atoms between the solitons. Because of exchange of 
atoms the motion of solitons in figures 3 is asymmetric in general except for 6 = tc and 
0. 

In figure 3 (a) for 6 = it in the interval of time 400, the solitons moved from y = ±15 
to ±38, respectively, corresponding to a final separation between the two solitons of 76. 
In figure 3 (b) for 5 = 37r/4 in the interval of time 400, the solitons moved from y = ±15 
to 35 and —41, respectively, again leading to a final separation of 76. Although the 
repulsion is same in both cases, an asymmetry has appeared in figure 3 (b) due to 
inelastic collision with exchange of atoms. In figure 3 (c) for 5 = 7r/2 in the interval of 
time 400, the solitons moved from y = ±15 to 38 and —26, respectively, corresponding to 
a final separation of 64. As 6 is reduced further the similar trend is maintained: reduction 
in repulsion due to inelasticity (exchange of atoms between the two solitons). In figure 
3 (d) we see that for 6 = tt/A the two solitons maintain their identity, exchange atoms 
and the overall interaction continues to be repulsive. During the interval of time 400, in 
figure 3 (d) the solitons moved from y = ±15 to 33 and —16, respectively, corresponding 
to a final separation of 49. For a further reduction of 6 from 7r/4 the inelasticity and 
asymmetry of interaction reduces. The mutual overall interaction becomes less and less 
repulsive which finally becomes attractive as 5 ^ 0. In figure 3 (e) we show the relative 
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Figure 3. One-dimensional probability t) vs. y and t for a train of two vortex 
solitons each with n — — 0.4, L — 1 and 5 = (a) tt, (b) 37r/4, (c) 7r/2 and (d) 7r/4, (e) 
7r/8 and (f) 0. 



motion of two solitons for 5 = vr/S. The two solitons have a tendency to combine to 
form a single soliton before they lead to two separate solitons. In time 400, solitons 
moved from y = ±15 to 22 and —11, respectively, corresponding to a separation of 33. 
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Figure 4. Ratio Nn/N vs. t for different initial phase 5 in Eq. (|3.1|l for a soliton-train 
of two L — 1 vortex solitions initially at y — i/q = ±15. Nfi is the number of atoms in 
the right soliton at = +15 at t = and the total number of atoms. 



The relative movement between the sohtons has reduced monotonically as 6 is changed 
from TT to 7r/8. From figure 3 (f) we see that for 6 = there is strong overall attraction 
and the two solitons interact and coalesce to form a single soliton which eventually gives 
birth to two solitons of same shape as the initial ones. The final separation between 
the two solitons in figures 3 (a), (b), (c), (d) and (e) for 6 = vr, 37r/4, 7r/2, 7r/4 and vr/S, 
respectively, are 76, 76, 64, 49 and 33 showing a gradual reduction in repulsion as S 
is reduced from vr to vr/S. In all cases the final separation is larger than the initial 
separation 2?/o(= 30) between solitons refiecting an overall repulsion. The asymmetry 
in the position of the two final solitions is zero for 6 = and vr and largest for S = n/A. 

Two features of figures 3 deserve some comments. First, the value of the phase 6 for 
which the interaction between two solitons is attractive. In the one-dimensional model 
study by Salasnich et al. it is attractive for all cos 5 > in agreement with a previous 
analysis [26] whereas in the present three-dimensional study with an additional radial 
trap there is no sharp transition from attraction to repulsion for a fixed 6. The presence 
of the radial trap in our study makes a slow transition from repulsion to attraction with 
the change of 6. However, there is clear attraction for 5^0 and repulsion for cos 5 < 
in both models. Further studies are needed for a complete understanding of the present 
three-dimensional solitons under transverse confinement. This transverse confinement 
is absent in the one-dimensional model of Salasnich et al. J2] • It is true that in the one- 
dimensional model of Salasnich et al. the radial trap frequency enters as a parameter [see 
their equation (1)], it does not include real three-dimensional nonlinear dynamics. The 
present three-dimensional solitons cannot be considered to be classic one- dimensional 
integrable text-book solitons of the type considered in references jl2|l26j. Hence the 
text-book wisdom is not directly applicable to the present analysis. 

The second feature is the asymmetry in figures 3 for an arbitrary 6 because of 
inelastic exchange of atoms between the two solitons which is absent in the one- 
dimensional model Although we report the results for angular momentum L = 1 in 
this study we have verified that similar asymmetry is present for L = 0. The asymmetry 
is also present in the analysis of ^14j of two interacting solitons. One of the two equal 
solitons increases and the other decreases in size due to inelastic exchange of atoms. 
In figures 3 the left soliton grows in size and moves slower, the right soliton becomes 
smaller in size and moves faster. The direction of particle exchange is determined by 
the way the phase is introduced between the two solitons in 1)3.11) . We verified that as 
6 is changed to —S the role of asymmetry is reversed. For a positive S due to particle 
exchange in figures 3 the left soliton grows in size. For a negative 6 the right soliton 
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grows in size in a symmetric fashion so that P{y,t, —6) = P{—y,t, 6). The asymmetry 
of interaction is absent in figures 2 and 3 (a) where the interaction is purely repulsive 
with no exchange of atoms. 

To study systematically the asymmetry of the interaction beween two solitons we 
plot in figure 4 the ratio Nn/N as a function of time for the cases presented in figures 
3. Here Nji is the number of atoms in the right soliton and N the total number in two 
solitons. The variation of Nji/N is quite similar to that found in ^3]. For S = n/A in 
both studies N^/N first attains a minimum and then oscillates, however remaining less 
than 0.5. A similar behavior is noted in the present study for 6 = tt/S. For S > 7r/4, 
Nji/N attains a minimum and then increases a little to attain a constant value smaller 
than 0.5. However, a quantitative comparison between the two studies is not to the point 
as the initial conditions of the two studies are different, e. g., the initial size and angular 
momentum L and the initial separation between the solitons. The details of interaction 
dynamics should depend on these initial conditions. In the present simulation we find 
that, for the change 6 —S, Nn/N Nl/N, where Nl = {N — Nji) is the number of 
atoms in the left soliton. 

Next we exhibit the profile of two interacting solitons undergoing particle exchange. 
From figure 4 we find that 6 = n/A leads to a large exchange of particles. For this purpose 
we consider the situation depicted in figure 3 (b) corresponding to two L = 1 vortex 
solitons of 72 = —0.4 each which are set at ?/ = ±15 at time t = with a phase difference 
of (5 = 7r/4 according to ()3.1|) . The profiles of the system at times t = 0, 100,200, 300 
and 400 are shown in figures 5 (a), (b), (c), (d) and (e), respectively. Due to exchange 
of particles the left soliton gradually becomes larger and larger consistent with figure 4. 
At t = 100 the solitons come closer, interact strongly by exchanging atoms and become 
wider and deformed. Then they gradually move apart and become narrower. Finally, 
at t = 400 they re-acquire Gaussian shapes. However, at t = 400 the right soliton is 
wider and shorter and accommodates a smaller number of atoms consistent with figures 
3 (b) and 4. At that stage the left soliton is narrower and taller and contains a larger 
fraction of the atoms. 

The interaction between more than two solitons with a general phase 6 among 
different pairs is much too complicated to be studied exhaustively here. If the phase 
difference 6 between two neighboring solitons is not close to zero, they experience overall 
repulsion and stay apart. However, for 6 close to zero they interact attractively and 
often a soliton could be lost as observed in the experiment of Strecker et al. 4J. This is 
illustrated in figure 6 where we consider three solitons all with same phase (5 = 0). Due 
to attraction the three solitons come closer, interact and form two solitons which after 
some time form a single soliton. This single soliton next decays to two, recombines to 
one, and decays to two again and never three original solitons are recovered. It is also 
possible that for some values of nonlinearity n and initial separation 2?/o two solitons 
with relative phase 5 = may coalesce to form a single soliton without ever decaying 
to two solitons again. These could explain a missing soliton observed in the experiment 
by Strecker et al. |4j. 

Throughout this investigation in the interaction of two equal solitons we assumed 
that the nonlinearity |?t.| for each is less than |?T,cr|/2, so that a stable solitonic condensate 
with total \n\ < |ncr| exists when the two coalesce. However, if two solitons each with 
\n\ > |ncr|/2 encounter for 6 = 0, the system is expected to coalesce, collapse and emit 
atoms via three-body recombination [S^. It is possible that in this case only a smaller 
single soliton survives. This might also explain some missing soliton(s) in experiment. 

In the investigation by Al Khawaja et al. [TT], to justify the repulsion between 
solitons in a train |3] a phase difference of 6 = tt was suggested between neighbors. 
They made a model for the evolution of S along the axial direction which they used for 
explaining the formation of a soliton train, with neighboring solitons of phase difference 
TT, after changing the scattering length from repulsive to attractive in a BEG as in 
the experiment i;4]. Such an order in phase seems to be not necessary for an overall 
repulsion between solitons. Almost any 6, except 6 near 0, is found to lead to overall 
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Figure 5. The tlic-dimensional wave function profile \ip{x,y)/x\ vs. x and y of two 
equal L = 1 rotating bright solitons of figure 3 (b) with n = —0.4 each and S = n/A at 
times (a) t^Q,{h)t = 100, (c) t = 200, (d) t = 300, and (e) t = 400. 



repulsion. In reference ^2], using an approximate analysis, it has been found that 
repulsive interactions require cos 5 < 0. In the present three-dimensional analysis, we 
find that attraction starts at 5 = 7r/2 and increases as 5 is reduced. However, the overall 
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Figure 6. One-dimensional probability P{y,t) vs. y and t for a train of three vortex 
solitons at and ±20 each with n = —0.4, L = 1 and 5 = 0. 



interaction remains repulsive until a small value of 6 close to zero is attained, when the 
solitons attract, lose their separate identity and coalesce. 

4. Conclusion 

To conclude, employing numerical solution of the GP equation with axial symmetry, we 
have performed a realistic mean-field study of interaction among two vortex solitons in a 
train and find the overall interaction to be repulsive except for phase 6 between neighbors 
close to 0. For phase 6 = tt between neighboring solitons in a train, they are found to 
repel and move away without exchanging atoms for both normal (L = 0) and vortex 
(L = 1) solitons. For phase 5 = in a two-soliton train, the solitons attract, collide, 
coalesce, and eventually come out. For other 6, overall repulsion prevails. However, 
there is an inelastic exchange of atoms between two solitons resulting in a change of size 
and shape. These unequal solitons travel in general with different speeds: the smaller 
soliton travels faster and the larger soliton travels slower. The unequal speed leads 
to an asymmetry in final positions of the solitons. By changing the sign of the phase 
between the two solitons the asymmetry in the final size and shape of the two solitons 
can be reversed. Except in the 5 ~ case, the solitons in a train stay apart and never 
cross each other as observed in the experiment by Strecker et al. [4j. For 5 ~ a single 
soliton can often disappear as a result of the attractive interaction among solitons, as 
observed experimentally by Strecker et al. jlj. These features of soliton interaction are 
present for both normal (L = 0) ^1] vortex (L = 1) solitons. Although the present 
study is performed in the absence of an axial trap as in one dimension 0, these general 
conclusions should remain valid for a weak axial trap as in the experiments jHIHl- The 
L = 1 vortex solitons can accommodate larger number of atoms and the present study 
may motivate future experiments with them. 
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